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A systematic study is performed by the ab-initio density functional theory of the anisotropy of the 
orbital moments (Iz) and the magnetic dipole term (Tz) in bulk Cr02. Two different band-structure 
techniques are used (FLAPW and LMTO-ASA), and the electronic correlations are treated by 
the local-spin-density approximation (LSDA), the LSDA-I- orbital polarization method, and the 
LSDA+(7 method. The calculated anisotropics of (Iz) and (Tz) are very large compared to Fe, Ni 
and Co but still a factor of 5 and 2 smaller than the anisotropics obtained from a recently suggested 
analysis of the X-ray magnetic circular dichroism spectra for a thick layer of Cr02. 
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Recent research within the fields of magnetic tunnel- 
ing and spin injection involves Cr02 as a promising ma- 
terial for electrodesi. Besides its potential importance 
for appUcations, this material also exhibits interesting 
physics originating from its half metallic nature and con- 
sequently the ferromagnetism due to the double-exchange 
coupling^, as well as from its orbital magnetism^ which 
is related to the spin-orbit coupling and the electronic 
orbital correlation effects. Consequently, there are al- 
ready several experimental^ii and theoretica P'^'^'^ inves- 
tigations of the orbital moments (Iz) in Cr02. 

However, much less information is available on the 
anisotropy of the orbital moments, i.e., on its dependence 
of the orientation of the sample magnetization (in the 
following the z axis of the external coordinate system is 
always chosen to be parallel to the magnetization direc- 
tion). It has been pointed out first by Bruno^ and later 
worked out in more detail by van der LaanS that the 
anisotropy of the orbital moment is closely related to the 
magnetocrystalline anisotropy energy. Furthermore, van 
der Laan has shown that there is an additional contribu- 
tion to the anisotropy energy arising from the anisotropy 
of the magnetic dipole term Tj , which is the expectation 
value of the magnetic dipole operator: 
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= -[cr-3r(f-cr)]. 



(1) 



In eq.(l) f is the unit vector in the direction of the posi- 
tion vector r, and cr is the vector of the Pauli matrices. 
In fact, it turns out (see below) that for Cr02 the two 
contributions are very large and of similar magnitudes 
but opposite in sign. 

A suitable method to investigate the anisotropy of (Iz) 
and (Tz) is the angle-resolved variant of the X-ray mag- 
netic circular dichroism^- (XMCD) which measures the 
XMCD spectra for various orientations of the sample 
magnetization. The orbital moments (Iz) then may be 
determined directly from the application of the XMCD 
orbital sum rule^^. In contrast, the application of the 
spin sum rul»l^ yields a combination of the spin moment 



(az) and of the {Tz) term. A separation of these two 
contributions is possible by applying the spin sum rule 
to XMCD spectra for various orientations of the sample 
magnetization. A precondition of the use of this method 
is that '^ij^{Tz)a ~ 0, where a represents one of the per- 
pendicular directions in which the z axis is oriented. It 
has been outlined in Ref. Q that the validity of this re- 
lation is not always guaranteed (for the case of Cr02 see 
below). 

Goering et alA applied the technique of angle-resolved 
XMCD to investigate the anisotropy of {Iz) and {Tz) for 
the Cr atom in a thick layer of Cr02 on a Ti02 substrate 
(which is required to stabilize the rutile structure). For 
Cr02 the analysis of the XMCD spectra of the L2 and L3 
edges of the Cr atom (corresponding to the 2pi/2 3d 
and 2p3/2 ^ 3d transitions) is more difficult because 
the spin-orbit coupling of the Cr 2p electrons is rela- 
tively weak so that the 2pi/2 and 2^3/2 levels are not 
well separatedi^. The subdivision of the XMCD spectra 
into the L2 and L3 contributions (which is required for an 
application of the spin sum rule) is therefore highly prob- 
lematic due to a possible quantum-mechanical mixture of 
the 2px/2 and 2p3/2 levels and/or due to a strong overlap 
of the two contributions on the energy scale. In their 
original analysis Goering et al.^ neglected the quantum- 
mechanical mixture. Furthermore, they separated the 
L2 and L3 contributions by an empirical (and theoretical 
not yet justified) extension of the van der Laan's method 
of the moment analysisi^ (which was proposed for sys- 
tems with very small crystal-field splitting of the valence 
states) to situations with large band splitting. Based 
on these two assumptions, they found extremely large 
anisotropics of {Iz) and {Tz). For instance, the differ- 
ence in the orbital moment for the magnetization along 
the c axis and along the a axis of the rutile structure 
was 0.083 fiB, which is extremely large compared to bulk 
materials with cubic symmetry (10~^ ^b, see Ref. IT^ 
and comparable to the anisotropy found for extremely 
thin Co layersiSi. The anisotropy of {Tz) appeared to be 
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about 0.2 ^B, two orders of magnitude larger than the 
absolute value of (T^) in Fe (0.004 ^b), Ni (-0.004 ^b) 
and Co (0.002 /iB)"®- For the true spin moment of the 
Cr atom a value of 1.2 /iB was obtained. To explain the 
magnetic moment of about 2 /iB per unit cell which was 
obtained by a SQUID measurement Goering et alA as- 
sumed a very large spin moment of about 0.4 /xb per O 
atom which they tried to explain in terms of a hybridiza- 
tion between chromium and oxygen. 

The question is whether the extreme anisotropics of 
{Iz) and {Tz) found by Goering et alA are indeed intrinsic 
properties of bulk Cr02 with rutile structure. Remember 
that their measurements were performed on thick layers 
of Cr02 on Ti02. So far there is only very little infor- 
mation on the structure of this layer, i.e., on possible 
oxygen deficiencies. Furthermore, it is not known how 
the two basic assumptions for the application of the spin 
sum rule (see above) affect the results for {Tz) and for 
the true spin moment (cr^). Indeed, most recently^S, an 
empirical method has been developed to take into ac- 
count approximately the effect of mixing of the pi/2 and 
P3/2 core states for the analysis of the XMCD spectra 
by the spin sum rule. The data of Goering et al. were 
reanalyzedSSi by this method, yielding a Cr moment that 
agreed reasonably well with the Cr moment suggested 
by the SQUID measurements after neglecting the O mo- 
ment, but the anisotropy of the (Tz) remained to be very 
large when obtained by the help of the above discussed 
moment analysis. 

In the present paper we investigate theoretically the 
anisotropy of (Iz) and (Tz) in bulk Cr02 with rutile struc- 
ture by calculations based on the ab-initio density func- 
tional electron theory. In the literature there are already 
some single results in this direction which - however - are 
in part contradictory^iSii. We therefore have performed 
a systematic study, which is based on two different band- 
structure techniques and which takes a special care of 
the electronic orbital correlation effects (again by two 
different methods) which are not included in the com- 
monly used local-spin-density approximation (LSDA) of 
the density functional theory. 

We performed comparative calculations by two dif- 
ferent band structure methods, the Wien97 code2i, 
which adopts the full-potential linearized-augmented- 
plane-wave method^Si (FLAPW) and the linear-mufhn- 
tin-orbital method in the atomic-sphere approximation^"^ 
(LMTO-ASA). Whereas in the first method the effec- 
tive crystal potential is treated exactly, it is spheri- 
cally averaged in each atomic sphere after each itera- 
tion step of the self-consistency cycle in the LMTO-ASA 
method. The exchange-correlation potential was calcu- 
lated in LSDA^. The electronic orbital correlation ef- 
fects, which are not included in LSDA, were taken into 
account by the OP termS^ or by means of the LDA-I-J7 
schemeSSiSi. For [/ = 0, the LT)A+U calculation is equiv- 
alent to a LSDA calculation. In the case of J7 7^ the 
exchange-interaction parameter J, which appears in the 
LDA+[/ scheme in addition to the parameter U, was 



fixedS to 0.87 eV, whereas U was an open parameter 
(the calculated^ screened value of U for Cr in Cr02 is 
U — 3 eV) . In the literature arguments for and against 
the use of the LDA-I-J7 method for the case of Cr02 are 
given (see, e.g., Ref. |5j and Refs. therein). Therefore we 
performed a comparative study based on various calcula- 
tional schemes for the electronic correlations. The spin- 
orbit coupling, OP term and hDA+U scheme were imple- 
mented in the LMTO code in Ref. |2^ The magnetization 
was set along the c or the a axis of the rutile structure. 
The calculations were performed for the experimental lat- 
tice parameters^- a — 0.4419 nm, c = 0.2912 nm and 
u = 0.303. 

The results are presented in Table 1. Hund's third rule 
is obeyed for both magnetization directions at any value 
of U because the orbital moments are always antiparal- 
Icl (parallel) to the spin moments at the Cr(0) sites, in 
agreement with the experimental observations^. 

We first discuss the orbital moments. From the 
LSDA calculation, the LSDA-f OP calculation and for the 
LDA+U calculation with U < 3eV we find Cr orbital mo- 
ments of about —0.05 /^b which magnitudes are smaller 
than the experimental values^ for Fe (0.086 /ib) and Co 
(0.153 ^b)- For the magnetization along the c axis our 
FLAPW calculation with [/ = 3 eV yields orbital mo- 
ments of —0.064 fiB and —0.0025 /iB at Cr and O sites, 
in a very good agreement with the XMCD results of Ref. 
i which give -(0.06 ± 0.02) /iB and -(0.003 ± 0.001) mb, 
respectively. We take this as a hint that U — 3 eV (which 
is also the calculated^ value for the screened U) is a good 
choice. This is further underpinned by the fact that our 
LSDA-I-OP calculation yields very similar results for the 
orbital moments of chromium and oxygen as the LDA-I-J7 
calculation for U between 2 and 3 eV. For U < 3 eV 
our FLAPW results for the Cr (O) orbital moment agree 
very well (perfectly) with those of a full-potential LMTO 
(FP-LMTO) calculation^. The probably less-accurate 
LMTO-ASA calculation yields larger Cr orbital moments 
than the FLAPW calculation. For f7 = 3 eV the LMTO 
calculation gives a value of —0.090 /iB for the c-axis ori- 
entation which agrees very well with the XMCD result of 
Ref. i4. For the magnetization along the a axis (for which 
there is no FP-LMTO result in Ref. 0) both our FLAPW 
and our LMTO-ASA calculations yield smaller orbital 
moments than those for the c-axis orientation, in agree- 
ment with the XMCD resultsi. The difference \ A{lz) \ = 

|(^z)c axis - (^z)a axisi IS about 0.015 /iB at J7 = 3 cV, 

which is indeed very large in comparison to the orbital- 
moment anisotropy of cubic bulk materials (10~*/iB), but 
still a factor of about 5 smaller than the orbital-moment 
anisotropy of Cr02 discussed by Goering et al^. In the 
FLAPW calculation the value of |A(Z2)| decreases dras- 
tically when reducing U and it amounts to 0.002 /tB for 
the LSDA calculation (the FLAPW calculation™ based on 
the generalized-gradient approximation (GGA) gives an 
orbital moment anisotropy of 0.001 /iB). In the LMTO- 
ASA calculation the reduction with decreasing U is less 
dramatic: for U = OeV we find |A(^2)| = 0.01 /iB (which 
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TABLE I: The calculated values of {Iz), {(^z), (Tz) and the anisotropies | A(Zz) | and | A(rz) |, aU in /ib, from the LSD A calculations 
{U — eV), the LDA-I-J7 calculations for various U and for the LSDA-I-OP calculations. The values in brackets are from the 
LMTO-ASA calculations, all other values are from the FLAPW calculations. 



is a factor of about 5 smaller than the value found by a 
LMTO-ASA calculation by Uspenskii et al.^). Increasing 
U beyond 3 /^b enhances the orbital moment anisotropy 
to the values comparable to 0.083 /ub discussed by Goer- 
ing et alA, but we doubt that such large values of U are 
physically reasonable for the Cr atom in Cr02. 

The magnetic dipole term (T^) at the Cr atom de- 
pends only slightly on the way we take into account the 
correlation-effects. The LSDA, LSDA-^OP and LDA-f 
([/ < 3 eV) calculations with the FLAPW method yield 
values of about —0.06 ^b (+0.04 ix-q) for the c-axis (a- 
axis) orientation. These values are large compared to 
the LSDA valuesi2, for Fe (0.004 /zb), Ni (-0.004 /xb) and 
Co (0.002 /xb), and they arise from a large crystal-field 
anisotropy. Such large values usually appear in systems 
with strongly reduced dimensionalitjsi^, and they show 
that the Tz-term cannot be neglected in the analysis of 
the XMCD spectra via the spin sum rule^^ when one 
wants to arrive at realistic values for the true spin mo- 
ment. The {Tz) values from the LMTO-ASA calculation 
are a bit smaller which results from the fact'^^ that in 
this method the spin and charge densities are calculated 
for an effective potential, which is spherically averaged 
in each atomic sphere. The anisotropy |A(Tz)| is very 
large, about 0.1 /iB, in the FLAPW calculation but still 
a factor of 3 smaller than the one discussed by Goering 
et alA. 

Finally, we have checked by the FLAPW method the 
applicability of the angle-resolved spin sum rule analy- 
sis for the case of Cr02. This analysis is appropriate if 
'^J2jTz}c./{az) < L For J7 = 3 eV (0 eV) we find a 
ratio of 0.09 (0.01)). 

Based on our values for |A^z| and lAr^l we calculated 
according to Ref. ^ the respective contributions to the 
magnetocrystalline anisotropy energy. They appear to 
be very large, of the same order of magnitude but op- 
posite in sign. Because small uncertainties in either of 



the two contributions induce very large errors for the 
magnetocrystalline anisotropy, it does not make sense to 
calculate the latter quantity from our data. 

Finally, we consider the results for the spin moments. 
In all our different types of calculations we find a Cr spin 
moment which is only very slightly anisotropic and which 
exhibits a magnitude of about 2 /iB, the value predicted 
by Hund's first rule for the case of a Cr^+ ion but con- 
siderably larger than the value of 1.2 /xb found by the 
special type of analysis of the XMCD data^ discussed 
above. For the O atom our FLAPW calculation yields 
for U = 3 eV a spin moment of about 0.1 /j,b which is a 
factor of 4 smaller than the one obtained by Goering et 
al^ with their special analysis of the XMCD results. Our 
FLAPW spin moments for oxygen are a bit larger than 
those from the FP-LMTO calculation of Ref. B Our 
LMTO-ASA data are even smaller than the FP-LMTO 
values. This probably results from the fact that sizes of 
the spheres over which the spin density is integrated are 
different in various calculations. The atomic sphere in 
the LMTO-ASA calculation is larger than the muffin-tin 
spheres in the FP-LMTO and FLAPW calculations and 
probably contains already a part of the Cr spin-density, 
which has an opposite sign. 

To conclude, we performed a systematic investigation 
by the FLAPW and the LMTO-ASA methods for the 
anisotropies of the orbital moment {1^) and the magnetic 
dipole term (T^) in Cr02. The electronic correlation ef- 
fects were taken into account by three different methods, 
the LSDA, LSDA+OP and the LDA+[/ method. The 
calculated anisotropics of (Z^) ( {Tz)) are very large but 
still a factor of 5 (2) smaller than those obtained by Go- 
ering et alA from their special type of analysis of XMCD 
data for thick layers of Cr02 on Ti02. To arrive at such 
large anisotropies we had to insert values for the param- 
eter U of the LDA-t-J7 method which are considerably 
larger than 3 eV and which seem to be physically unrea- 
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sonable. The only chance to arrive at larger anisotropics 
might be to abandon the mean-field approximation which 
is adopted in the LDA+C/ method for the electronic on- 
site correlations and to perform a calculation within the 
dynamical mean field theory (DMFT)'^^. 

On the experimental side, more investigations are re- 
quired on the real structure of the considered layers as 



well as on the way how to analyze the XMCD spectra 
which arise from energetically very close core levels. 
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